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Abstract 

If there was a first order phase transition in the early universe, there should be an 
associated stochastic background of gravitational waves. In this paper, we point out 
that the characteristic frequency of the spectrum due to phase transitions which took 
place in the temperature range 100 GeV - 10'' GeV is precisely in the window that 
will be probed by the second generation of space-based interferometers such as the Big 
Bang Observer (BBO). Taking into account the astrophysical foreground, we determine 
the type of phase transitions which could be detected either at LISA, LIGO or BBO, in 
terms of the amount of supercooling and the duration of the phase transition that are 
needed. Those two quantities can be calculated for any given effective scalar potential 
describing the phase transition. In particular, the new models of electroweak symmetry 
breaking which have been proposed in the last few years typically have a different Higgs 
potential from the Standard Model. They could lead to a gravitational wave signature 
in the milli-Hertz frequency, which is precisely the peak sensitivity of LISA. We also 
show that the signal coming from phase transitions taking place at T~ 1-100 TeV could 
entirely screen the relic gravitational wave signal expected from standard inflationary 
models. 



1 Introduction 



Direct detection of gravitational radiation will hopefully soon become a reality, with the 
operation of the first generation of interferometers such as (kilometer-scale, ground-based) 
LIGO [1] and VIRGO [2] and (million- kilometer-scale, space-based) LISA [3]. Those in- 
struments will allow us to probe gravitational waves produced by astrophysical objects at 
relatively low red-shift (black hole binaries, neutron star binaries, white dwarf binaries, su- 
pernovae, pulsars...) In addition, gravitational waves (GW) can provide information about 
particle physics at unexplored high energies. The weakness of the interaction with matter 
is a major obstacle for detection of gravitational waves but it also has the virtue that the 
information they carry about the state of the universe at the moment of their production 
has been unaltered. They are precious information on the very mechanism that produced 
them. GW can be produced by core collapse of supernovae, first-order phase transitions, 
vibration of cosmic strings, preheating, dynamics of extra dimensions. . . 

Among those well-motivated but hypothetical cosmological sources of GW, there is at 
least one that we are convinced exists: the GW produced during inflation. This signal is 
expected to be very tiny. Quantum fluctuations in the inflaton field during infiation leaves 
behind a residue in density perturbations observed in the Cosmic Microwave Background 
(CMB). They also lead to a background of GW whose properties couple with those of 
density fluctuations. As the CMB anisotropics are affected by GW, the WMAP constraint 
on the energy scale of inffation fixes a bound on the size of the GW signal due to infiation 
^Gwh^ ^ 10^^"'"^ — 10^^'' [4]. This is several orders of magnitude below the best sensitivity 
of the first-generation of interferometers. However, attempts to detect this relic primordial 
background are very strongly motivated. This is a main goal of the second generation of 
space interferometers, in particular, the Big Bang Observer (BBO) [5], the follow-on mission 
to LISA, which would become a reality within twenty or thirty years (by comparison, LISA, 
if funded, should be operational by 2014). 

The present work focusses instead on the detectability of GW from first-order phase 
transitions. The corresponding relic GW background encodes useful information on these 
major symmetry-breaking events which took place in the early universe. In contrast with 
the infiationary spectrum, the spectrum is not fiat, with a characteristic peak related to 
the temperature at which the phase transition (PT) took place. This signal can actually 
be higher by several orders of magnitude than the signal expected from infiation and in 
some cases can entirely screen it. One symmetry- breaking event which for sure took place 
in the early universe is electroweak (EW) symmetry breaking. What we do not know yet is 
whether it was a first order phase transition, in which case it proceeded through nucleation 
of bubbles resulting in a large departure from thermal equilibrium. Bubble collision and 
associated motions in the primordial plasma are sources of gravitational waves. The charac- 
teristic frequency of the signal is close to the Hubble frequency at the time of the transition 
H{Tew) ~ lO"-*^^ GeV. Once redshifted to today, this corresponds to mHz frequencies, which 
is precisely the frequency band that LISA is sensitive to. It is therefore very exciting that 
LISA could help providing information on the EW scale, in particular on the nature of the 
EWPT. 

The GW spectrum resulting from first order PT was computed in the early nineties [6-9] 
but this topic has not received much subsequent attention, as it was found out that there 
is no first order EWPT in the Standard Model given the experimental bound on the Higgs 
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mass [10]. It was realized ten years after the original calculation of [6-9] that turbulence 
in the plasma could be a significant source of GW in addition to bubble collisions [11,12]. 
Subsequently, the authors of [13] studied the GW signal due to a first order EWPT in 
the Minimal Supersymmetric Standard Model (MSSM) and its NMSSM extension. Finally, 
Nicohs [14] did a model-independent analysis for the detectability of GW with LISA. 

We believe that it is time to revisit this question for two reasons: The nature of the 
EWPT will start to be probed experimentally at the LHC. Indeed, it depends essentially on 
the Higgs sector of the theory or any alternative dynamics for EW symmetry breaking. In 
the last few years, new models of EW symmetry breaking have been suggested (little higgs, 
gauge-higgs unification, composite higgs, higgsless models ...) and the nature (smooth cross- 
over or first-order) of the EWPT in these new frameworks remains unknown. Second, the 
technology for gravitational wave detectors has made advances [15] and we think it is timely 
to redo a model-independent analysis not only for LISA but also other devices. 

LIGO is sensitive to much higher frequencies (from a few Hz to a few hundreds of Hz) thus 
it is in principle sensitive to phase transitions which took place at much earlier epochs. For 
instance, we will show that if there was a very strong first order PT at temperatures of order 
10^ GeV, the ultimate stage of LIGO (LIGO-III, correlated) could detect the corresponding 
peak (as already pointed out in [6]). The second generation of interferometers will be able to 
say much about the possible existence of early universe first-order phase transitions. Indeed, 
BBO will be sensitive to signals from PT which would have taken place in the temperature 
range T ~ 100 GeV- 10^ GeV, even if not necessarily exceptionnally strong. 

In this paper, we start with some generalities on stochastic GW backgrounds including 
the astrophysical background. We also recap what would be the observable redshifted signal 
we would observe after the GW have propagated forward from the phase transition until 
today. Section 3 reviews the key formulae used in the theoretical predictions of the GW 
spectrum due to first-order phase transitions. There is nothing new in this part. However, 
this formalism had so far only been exploited to study the detectability at LISA of GW due 
to a first-order electroweak phase transition. In section 4, we apply it to any other phase 
transitions taking place at higher temperatures and compare them with the sensitivities of 
not only LISA but also LIGO and BBO. Predictions can be presented in a model-independent 
way as a function of two quantities, namely a (~ latent heat) and (~ duration of the 
phase transition), which can be computed for any given effective scalar potential describ- 
ing the transition. For each temperature, we identify which values of a and f3 lead to an 
observable signal. Particle model builders can then test their favourite scalar potential by 
computing its corresponding values of a and (3 and see whether it can give rise to a detectable 
GW signal. We comment on some specific examples of particle physics models. 



2 Astrophysical versus cosmological GW background 

Stochastic backgrounds are random gravitational waves arising from the incoherent super- 
position of a large number of independent, uncorrelated sources that cannot be resolved 
individually. They are discussed in terms of their contribution to the universe's energy 
density, over some frequency band: 

naw(f) = —P^ (1) 
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By their very nature, stochastic GW are indistinguishable from the detector noise. Ground- 
based detectors look for them by coordinated measurements (comparing outputs of multiple 
detectors to find sources of correlated noise) while LISA can extract the instrumental noise 
power by combining the signals from its three spacecrafts. For technical aspects related to 
the detection of a gravitational wave stochastic background, see Ref . [16] . 

2.1 Astrophysical foreground 

Searching for GW waves of cosmological origin is an ambitious goal. There is a huge fore- 
ground due to astrophysical sources which in principle makes detection impractical. Once 
the signals from every merging neutron star and stellar mass black holes have been identified 
and substracted, the primary sources of foreground signals are galactic and extragalactic 
binaries. The galactic background produced by binary stars in the Milky Way is many 
times larger in amplitude than both the extragalactic foreground and LISA's design sensi- 
tivity. However, it can be substracted because of its anisotropy, being mostly concentrated 
in the galactic plane. Irreducible background comes from extragalactic binary stars and is 
dominated by emission from white dwarves (WD) pairs. The corresponding GW spectrum 
was estimated in Ref. [17] where limits are placed on the minimum and maximum expected 
background signals. Ref. [17] points out that at frequencies / ^ 50 mHz, there will be too 
many individual WD-WD sources contributing in each resolution element to be completely 
resolved and substracted source by source by missions with plausible lifetimes. However, 
much of the fiux comes from relatively nearby sources, and the WD-WD numbers drop 
rapidly above 50 mHz. Thus it may be possible for future missions more sensitive than LISA 
to substract this background at high frequencies [17]. In our figures, we plot this background 
coming from unresolved compact white dwarf binaries assuming that it can be removed at 
frequencies above 50 mHz. At higher frequencies, the dominant foreground GW sources 
are inspiralling neutron star-neutron star, neutron star-black hole and black hole-black hole 
binaries. These have to be individually identified and substracted. This problem is discussed 
in [18] and in our BBO detectability analysis we optimistically assume that this foreground 
can be substracted. 

2.2 Relic background from cosmological processes 

The GW background due to early universe events is stochastic as the signal comes from 
the superposition of incoherent sources originating from a huge number of different horizon 
volumes. For instance, the size of the horizon at the time of the electroweak phase transition 
was much smaller than today (10~^^ GeV) ^, corresponding to a tiny fraction of degree on 
the sky today. Even if we assumed that there were two bubbles per horizon volume (actually 
there would typically be several hundreds of them as we will see later) , we would be unable 
to resolve the signal coming from their collision. The signal comes from bubble collisions 
which took place in many independent universes. 

We work in a standard Friedman-Robertson-Walker (FRW) cosmology, a{t) is the cos- 
mological scale factor. At the energy scales considered, we assume a radiation-dominated 
era. Gravity waves produced at with a characteristic frequency propagate until today 
without interacting. Their energy density redshifts as and their frequency as a~^. The 
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characteristic frequency we observe today is 

(n \^/^ T 

(2) 

where we used the adiabaticity of the expansion of the universe (meaning that the entropy 
per comoving volume S oc a?gs{T)T^ remains constant). 

9.(T)- E 9.i^i^l E 9.{f) (3) 

i=bosons i=jermions 

Qi counts the internal degrees of freedom of the i-th particle and the sum is over relativistic 
species. Today, 5's(To) ~ 3.91 (assuming three neutrino species) and Tq = 2.725/^ = 2.348 x 
10~^^ GeV. It is convenient to express the frequency in terms of the Hubble frequency at the 
time of GW production: 

/ ^ 6 X 10-=^mHz ( ^* A (4) 
VlOO; 100 GeV ^ ' 

The remarkable fact is that for T ~ 100 GeV and f*/H^ ~ 10^, (as expected for weak scale 
processes as will be explained below), the peak frequency of the GW spectrum is in the 
milliHertz, just in the band of LISA. 

The fraction of the critical energy density in gravity waves today is 

^cw = '-^ = new. - 1.67 X 10-/.- f now. (5) 

Pc \aoJ KHqJ \ 9* J 

where we used 

PGW = PGW* f— V , Pc = Pc*^ and Ho = 2.1332 xhx 10"^^ geV (6) 
\aoJ Hi 

is the number of relativistic degrees of freedom at which enters the definition of the 
energy density and not the entropy (it is given by Eq. Q where the cubic power is replaced 
by a quartic power), flew* is the fraction of energy density of the universe at the time of the 
transition which is in gravitational waves. The peak sensitivity of LISA would correspond to 
detect flawh"^ ~ 10~^^. This means that to detect a signal at LISA, we need flaw* ^ 10^^ 
while at BBO, we can probe smaller fractions, ^gw* ~ 10~^^ — 10~^. 

The remaining task is to estimate and fi*. Theoretical predictions of relic GW back- 
grounds are subject to large uncertainties which depend on the cosmological mechanism. 
However, we can get a reasonable estimate of the characteristic frequency, the form of the 
spectrum and the typical intensity. We will now review the main results in the case of GW 
produced during first order phase transitions. 



3 GW from first-order phase transitions 

Phase transitions are commonly described by the effective potential of the scalar field (ei- 
ther elementary or composite) responsible for the dynamics. First-order phase transitions 
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are triggered if there exists a temperature at which a barrier separates two degenerate min- 
ima. This happens for instance if there are negative cubic or quartic self couphngs for the 
scalar field. In this case, phase transitions proceed via nucleation of bubbles of the low- 
temperature phase within the high-temperature phase. Bubble nucleation occurs through 
quantum tunneling and thermal fluctuations. These bubbles then expand and merge, leav- 
ing the universe in the low-temperature phase (commonly the broken-symmetry phase). As 
a bubble expands, part of the liberated latent heat raises the plasma temperature while 
the other part is converted into kinetic energy of the bubble wall and bulk motions of the 
fluid. Because of its spherical symmetry, a single expanding bubble produces no gravity 
waves. Only after bubble collisions destroy the spherical symmetry is gravitational radia- 
tion emitted. High velocities and large energy densities provide the necessary conditions for 
producing gravitational radiation. There are two sources of gravitational waves: the actual 
collision of bubbles and the turbulence in the plasma due to bubble motion. The resulting 
spectrum of gravitational waves has been studied in details in [6-9, 11-14]. The turbulence 
spectrum was recently revisited in Ref. [19]. Re-examination of the bubble collision spectrum 
is underway [20]. 

Remarquably, these predictions only depend on the grossest features of the bubble col- 
lisions. Gravitational radiation is insensitive to the internal structure of colliding bubbles, 
in other words, to the small scale conflguration of the scalar held in the colliding region [6]. 
As conflrmed by numerical simulations [7] , the enveloppe approximation works very well. It 
consists in neglecting the dynamics of the collision (overlapping) region. Kinetic energy is 
supposed to be concentrated in the uncollided (but spherically asymmetric) bubble walls. 
Once bubble walls collide, they stir up the plasma at a scale comparable with their radii at 
the collision time, leading to turbulence which also induces gravitational emission. 

3.1 Key parameters characterizing the GW spectrum 

A crucial parameter for the calculation of the gravitational wave spectrum is the rate of 
variation of the bubble nucleation rate, called /3. This quantity flxes the characteristic scale 
in the problem, the size of bubbles at the time of the collision, and therefore the characteristic 
frequency The duration of the phase transition is given by f3~^ and the size of bubbles 
is typically Rb ~ Vhf3~^ where Vb is the velocity of the bubble wall. The initial size of the 
bubble at the time of nucleation (of the order of T~^) is negligible compared to (3~^ which 
is of the order of the horizon size. 

The second crucial parameter characterizing the spectrum of gravitational waves is a = 
e/ Prodi the ratio of the latent heat liberated at the phase transition (e=latent heat) to the 
energy density in the high energy phase, commonly being radiation energy density, e is not 
necessarily vacuum energy, see for instance [21]. a and j3 are evaluated at the nucleation 
temperature and determine entirely the GW spectrum. They can be computed once we know 
the effective action for nucleating bubbles ("critical bubbles") which can be computed for 
any scalar potential describing the phase transition. Therefore, given a scalar potential at 
finite temperature, K(0, T), this is enough to derive the predictions for the GW spectrum. 

The rate of bubble nucleation is 

V{t) = A(i)e-^W (7) 
The prefactor A{t) has units of energy to the fourth power, A{t) ~ A1^, where Al ~ T" is the 
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typical energy scale of the transition and S{t) ^ S-^/T is the euclidean action of the critical 
bubble. 



S's = / 4:Tcr dr 



is the free energy of a critical bubble, (pb is the bubble profile of the critical bubble obtained 
by solving for the "bounce" : 



dr"^ r dr 



dV 



0, with 



dr 



and 







(9) 



r=0 



Most of the time variation of V{t) is in S'(t), and [3 is defined as: 



dS_ 
It 



1 dT 
f~dt 



(10) 



where is the time when the transition completes. In a neighbourhood of t*, r{t) grows 
exponentially with time as S{t) = S{t^) — (3{t — t^) + .... From the adiabaticity of the 
expansion of the universe, dT/dt = —TH, where H is the expansion rate of the universe, 
and we obtain 



JL = T — 
* dT 



* dT [tJ 



fill 



The temperature of the transition T^, is defined as the temperature at which the probability 
for nucleating one bubble per horizon volume per horizon time approaches 1. This guarantees 
that bubbles percolate even if the universe is infiating. This translates into 



S 



-4 In 



T\_ 

nipi 



[12) 



f3/H is dimensionless and mainly depends on the shape of the potential at the time of 
nucleation. According to (fT^ . it depends only logarithmically on the energy scale. For a 
phase transition at the weak scale ~ 10^ and this justifies what we said after Eq. (jH) 
(where stands for While the size of the bubble increases by orders of magnitude 
between nucleation and percolation (its initial radius R ~ (p/VAV at the time of nucleation 
can be neglected; all what matters is the typical size at the end of the transition which is 
given by Rb ~ VbAt ~ VbP~^), (and thus a and /3/H) is essentially unchanged between 
nucleation and percolation. 

Note that the ratio f^/H also fixes the number of bubbles. The number density of bubbles 
is roughly f3~^r{t) so that the number of bubbles in one horizon volume is (3~^r/H^, which, 
at the end of the transition is of order H/ (5. 

The latent heat is the sum of two contributions. The first one is the difference in free 
energies between the stable and metastable minima (which vanishes at Tc) while the second 
one comes from the entropy variation As (which is non zero at Tc in a first-order PT). This 
leads to the following formula for e: 



-AF - TAs = {-^V + TdV/dT)r 



(13) 



To conclude, for the analysis of gravity wave emission, the only two relevant quantities 
are the amount of latent heat injected into the plasma, Eq. (fT^. and the nucleation rate, 
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Eq. (|7j). So once the bubble action is computed (search for the bounce solution Eq. 0), 
everything is known. The fact that the phase transition is described or not by a fundamental 
scalar field is irrelevant. 

In practise, instead of solving Eq. Q, one can use either the thin wall or thick wall 
limits to approximate the bounce solution. In the regime of large supercooling (which is the 
one of interest as far as large signals of gravitational waves are concerned), the thick wall 
approximation is adequate. In addition, if the temperature is decreasing and the ratio T/T^ 
is getting low (T^ is the critical temperature at which the free energies in the two phases are 
equal) one can use the 4D euclidean action 5*4 to evaluate the nucleation rate rather than 



3.2 Scaling expectations 



The energy density in gravitational waves coming from bubble collision can be estimated by 
naive dimensional analysis as follows. The quadrupole formula for the power of gravitational 
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emission is Pqw = ^HQij ) ) where G is the Newton's constant and Qfj is the quadrupole 
moment of the source which is T^^^, the transverse traceless piece of the stress tensor. We 
can write 



■■■TT mass of system in motion x (size of system)^ 

CJ ■ ~ 

^-^ (time scale of system)^ 



kinetic energy 
time scale of system 



(14) 



thus, Pgw ~ G^kin- Let n be the efficiency factor which quantifies the fraction of the 
vacuum energy which goes into kinetic energy of bulk motions of the fiuid (as opposed to 
heating) : 

Ekin K a Prad {VbP'^f (15) 



Using G ~ H^/ptot* and pt. 



ot* 



[1 + a) Prad we get {d/dt leads to a factor) pcw* 



Ecw/iviP 2) where Eqw = PgwP ^- Finally: 

Pgw* 



GW* 



Ptot* 




(16) 



To get a large signal, we need to be small and a to be large, in other words, the phase 

transition should last as long as possible and the latent heat should be maximized. The 
scaling obtained is very close to what a more rigorous calculation gives for the redshifted 
value of the energy density evaluated at the peak frequency [9]: 

^^coll \Jc( 



1.1 X lO'^K^ 



a 



1 + a 



vl 



0.24 + vl 



100 



g* 



1/3 



(17) 



while the peak frequency is [9] 

^ 5.2 X lO-^mHz 



JL 



lOOGeV 



100 



1/6 



(1^ 



The scaling is slightly different in the case of GW from turbulence in the plasma (the analysis 
of turbulent motions of [11] was generalized in [12,14]): 

n 2 r. ^^-1 1/3 



a..b /i'(^.b) ^ 1.4 X iQ-%ivi 



100 

g* 



(19) 



with [12, 14] 
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- 3.4 X lO-^mHz^ 

Vb 



JL 



lOOGeV 



100 



1/6 



(20) 



The turbulent fluid velocities Ug are smaller than the bubble expansion velocities, unless 
turbulence is extremely strong. As a result, according to (fT^ and (pH) . the peak for the 
turbulence spectrum will be shifted to lower frequencies. Note that the peak frequency of 
the collision signal does not depend on a. Because /t^.t ~ /con x (ug/vb) and {ug/vb) is an 
increasing function of a, as a increases, the collision peak gets hidden by the high-frequency 
tail of the turbulence signal^ . 

The signal from turbulence is more promising than the signal from bubble collision due 
to the different scaling with a. Indeed, the scaling with a of the velocities and the efficiency 
factor are : 



Vb[a) 



1/V3 + Ja^ + 2a/3 



1 + a 



[22] 



u.Aa) 



+ na 



[14] 



(21) 



K(a) 



1 + 0.715« 



0.715a + 




[9] 



(22) 



The formulae (fTTj) and ()19|1 for VLh'^ show that the amplitude of the signal does not 
depend on the energy scale of the transition but only on the dimensionless parameters a and 
H^:/i3, in other words on the shape of the scalar potential at the time of nucleation. Very 
roughly, taking P/H^ ~ Ss(T^)/T^ ~ (9(100), it is clear that we need a ~ C(l) if we want 
Qh^ > 10~^° (to see something at LISA). Our experience with possible values of a and 
is based on the studies of potentials describing the electroweak phase transition. The values 
of a and j3/H^, are related to the ratio 0(T^,) /T^:, another quantity characterizing the strength 
of the phase transition, where is the vacuum expectation value of the Higgs. For instance, 
in the Minimal Supersymmetric Standard Model (MSSM), a is typically smaller than 0.1 
while p/H^ is larger than 1000 [13]. On the other hand, in the NMSSM, the authors of [13] 
found values of a ~ 0{1). Values of a larger than 1 correspond to a phase transition which 
is so strong that it is at the borderline of not being able to take place because the free energy 
of a critical bubble is too large: either the barrier separating the two minima is too large 
(the barrier even exists at zero temperature) or the distance in field space between the two 
minima of the potential is too big. The same applies if f3/H^ is smaller than (9(100). This 
situation is encountered not only in the NMSSM but also in effective theories with large 
negative quartic couplings [23], as will be presented in [25]. Note that for a given scalar 
potential, a and are actually correlated as a large value of a will be associated with a 

small value of Indeed, a is proportional to the latent heat, which grows as AV^, the 

depth of the potential at the minimum, is getting bigger. On the other hand, S3/T (and 
thus P/H^) typically scales like 1/ \/ . 

Let us make a few comments concerning the situation where a ~ 1. In usual phase 
transitions, this means that the vacuum energy is of the same order as the radiation energy 
density, therefore inflation starts before bubbles percolate. However, if Eq. p2|) is satisfled, 
this guarantees that the phase transition can still complete. The number of e-foldings of 



^Although the relation /tm-b ~ /coii x {ug/vb) is likely to be revised in Ref. [20]. 
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inflation is \n{V^^^/T^), which is typically less than 1. However, if the transition is very 
slow, i.e. (3/H < 1, one should take formulae (|T7jl and (fTIHl with caution as they are derived 
neglecting the expansion of the universe, which is not a good approximation ii f3/H < 1. 

Fig. ^ shows some GW spectra illustrating the predictions for various temperatures with 
representative values of a = 0.4, 1 and jd / H^. = 100,800,3000. Those plots also exhibit two 
examples of signals from inflation (and taken from [4]) for comparison. The scale of inflation 
is constrained by the CMB to be Ej < 3.4 x 10^^ GeV [4]. This fixes the largest signal we 
could expect from inflation as flQ^h"^ oc Ej. We also include the signal corresponding to 
Ej = 5 X 10^^ GeV which could be observed at BBO. 

4 Scanning the (a, plane 

In our analysis below, we will use the formulae of the previous section as well as the fact 
that the spectrum flcouh"^ is expected to increase as /^'^ while at high frequencies it drops off 
as f~^'^ and that fiturb^^ increases as while at high frequencies it drops off as /~''/^. This 
is already all summarized in the letter [14]. However, Ref. [14] focuses on the detectability 
at LISA of GW from a T = 100 GeV phase transition. We are now using this formalism 
to look in more details at the detectability of GW coming from any other 1st order phase 
transitions at future interferometers. We repeat that we are working at the level of an- 
order-of-magnitude estimate. Magnetohydrodynamical effects could make the slope of the 
turbulence high frequency tail smaller [14] and in any case, for a more precise analysis, the 
calculation of the power spectrum should be revisited first. We compare the GW spectra 
resulting from PT occurring at temperatures in the range [100 GeV, 100 PeV] with the 
sensitivities of LISA, BBO and LIGO correlated third generation (and taken from [24]). The 
BBO sensitivity is approximate and may change in the final design. We are actually using 
the sensitivity of BBO Corr, its correlated extension, which correlates two detectors, namely 
two LISA-like constellations (each LISA-like constellation orbits around the Sun at 1 AU and 
consists in three spacecrafts in a triangular configuration) that will allow to do correlations 
to measure the stochastic background. As discussed in Section EIH we take into account as 
well the irreducible background due to extragalactic white dwarf (WD) binaries. 

For each temperature, we are making a full scan of the (a, /3/if*) parameter space and 
determine the regions where at least one of the peaks is observable. According to Eq. (|17| 
I18I19|20|) . various situations can arise: 

• For relatively low a, the turbulence and collision peaks are well separated and can be 
observed. This is the ideal situation as the observability of these two peaks would be 
a smoking gun for the phase transition origin of these GW. The ratio of the two peak 
frequencies is a predicted function of a. In somes cases, the turbulence peak is at too 
low frequency to be observed by LISA or BBO but the minimum separating the two 
peaks is visible. 

• At larger a (> 0.64), the collision peak is hidden by the high frequency tail of the 
turbulence peak. However, there is a characteristic change of slope in the high fre- 
quency tail. Depending on the temperature of the transition, this change of slope can 
be observed or not. 
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I}Gwh2 T=100GeV I}awh2 T=1TeV 




Figure 1: Spectrum of gravitational waves expected from a first order phase transition (solid 
blue line) for four temperatures and for some choices of {a, (3 / H) values. The dashed red lines 
are the (approximate) predicted sensitivities of LISA, BBO, LIGO-III. The horizontal dashed 
green lines are the gravitational spectra expected from inflation, for two scales of inflation, 
for comparison. The black dashed curve is the estimate for the irreducible foreground due 
to white dwarf binaries (from [17]). At large a, only the peak from turbulence can be seen 
as well as a change of slope (shown as a circled cross) corresponding to the high frequency 
tail of the bubble collision spectrum. For low a, it is possible to see the collision peak as 
well. 



Our contour plots show the region where the turbulence peak is observable and the region 
where either the collision peak or the slope change is visible, at LISA (Fig. BBO (Fig. E]) 
and LIGO (Fig. The vertical line separates the low a region where the two peaks are 
well separated from the large a region where only the change of slope is visible. The lower 
horizontal bound is due to the fact that we cut the sensitivity of both LISA and BBO at 
10~^ Hz. Because of this frequency limit in the sensitivity, we cannot probe phase transitions 
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Figure 2: Different configurations of the signal versus the instrument sensitivity to show the 
quahtative dependence on parameters. The upper blue region is where the turbulence peak 
is observable while the lower red one is the region where either the collision peak or the 
point of slope change is visible. Precise locations of these different regions depend on the 
experiment and the temperature of the transition as illustrated in Figs. El ID and El 



below a ~ GeV and thus the QCD phase transition. In the BBO plots, we show the very 
important effect of the WD foreground on the detectability at BBO. 

Note that we could have also made contours corresponding to cases where none of the 
peaks are observable but the high or low frequency tails can still be detected. This will 
clearly enlarge the detectability region and this is work in progress. 



LISA will be able to detect the peak of GW from a 100 GeV first order PT only if it is 
extremely strong (a > 0.5 or (3 / H < 1000). The Higgs potential of the MSSM does not 
satisfy this requirement but it can in the NMSSN [13]. Higgs potentials with negative quartic 
couplings can also trigger strong EWPT as was shown in [23]. The corresponding prospects 
for GW detection will be presented in details elsewhere [25]. There are also exciting large 
signals expected from the high temperature behaviour of a warped extra dimension [21]. If T 
is ~ 500 GeV rather than 100 GeV, the GW peak will coincide with LISA's best sensitivity 
frequency and a larger region will be detectable. The prospects for detection of GW from 
the 100 GeV EWPT are very good at BBO. At these low frequencies, the WD foreground 
is lower. For instance, the turbulence peak for a > 0.3, f3/H ~ 200, is above the WD 
foreground and can be seen by BBO. Note also that at T = 500 GeV, values of as large 
as 10^-10^ can be probed. 
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Figure 3: Contours delimiting the region in the {a, (3/H) plane for which there is an ob- 
servable peak at LISA. The upper blue region is for the turbulence peak while the lower red 
one is the region where either the collision peak or the point of slope change is visible. Left 
of the vertical green line, the collision peak is visible. 
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Figure 4: Same as Fig. IHl but for BBO. The effect of including the constraint from the 
irreducible WD foreground is displayed and limits the observable regions from the uncolored 
ones to the ones in plain colors. As the temperature increases, the peaks are shifted to higher 
frequencies, thus the effect of the WD foreground becomes less significant. 
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Figure 5: Same as Fig. 01 but for LIGO-III. 
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Figure 6: Below each line (each one associated with the temperature of the phase transition), 
the gravitational wave signal at BBO from 1st order phase transitions entirely masks the 
signal expected from inflation. This plot strongly depends on the scale of inflation, which 
was chosen here to be Ej = 3.4 x 10^^ GeV. 



4.2 T=l TeV 

It is quite exciting that a 1 TeV PT with a > 0.4, P/H ~ 200 can be seen by LISA. 
Recent radical proposals to address the hierarchy problem predict rich new phenomena at 
the TeV scale. For instance, new dimensions at a TeV could give rise to observable signals at 
LISA [21]. There was also a recent study of the high temperature behavious of Little Higgs 
Theories where it was shown that EW symmetry was restored precisely at a temperature 
of order ~ 1 TeV [26]. This transition appeared to be first order. A more detailed analysis 
would be required to determine whether this PT could be strong enough to lead to an 
observable spectrum of GW at LISA. Unfortunately, this takes place in the regime where 
the effective theory ceases to be under control. This is nevertheless an interesting prospect. 
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If 13 /H > 100, the turbulence peak is above the WD foreground for a > 0.2, and thus 
can be seen by BBO. In addition, BBO can see the high frequency tail of the collision peak 
which covers good part of, if not entirely, the inflation signal. 

4.3 T=10 TeV 

At these temperatures, the peak cannot be probed by LISA, unless (3 / H < 10^. On the 
other hand, LISA can still probe the low frequency tail of these spectra and is therefore a 
compelling tool to probe scales that LHC will not be able to reach. As the temperature 
increases, the peaks are shifted to higher frequencies, thus the effect of the WD foreground 
becomes less significant and quite weak first order phase transitions can be probed. And the 
high frequency tail of the collision peak can entirely screen the inflation signal, depending on 
the scale of inflation (see Fig.lHl). If ~ 1000, it is possible to see both the turbulence and 
the collision peaks for 0.1 < a < 0.64 and assuming that the inflationary scale is sufficiently 
low. For (3/H 200, the collision peak can be seen for a as low as ~0.05 if the inflationary 
signal is below the BBO sensitivity. 

4.4 T=100 TeV 

If (3 / H ~ 200 and a > 0.2 the inflation signal is for sure entirely covered. If the inflation 
scale is below ~ 5 x 10^^ GeV, a as low as 0.05 could be detected at BBO. Two peaks can 
be seen if 0.1 < a < 0.64 and P/H 200. At larger P/H only the turbulence peak will be 
seen. 

4.5 T=10^ GeV 

This is a particularly interesting the same signal could be observed by both BBO 

and Ligo-III. Specifically, a phase transition with a = 0.8 and ~ 200 would give a 

turbulence peak observable by Ligo-III while the low frequency tail would be observable by 
BBO. In this example, the infiation signal would be hidden except in a very narrow frequency 
range between 50 mHz and 80 mHz. 

Interesting signatures end at this energy scale. Phase transitions at T> 10^ GeV cannot 
be probed by any of the planned interferometers. 

5 Conclusion 

We have shown that the GW background from early universe phase transitions may become 
relevant for a second generation detector such as the Big Bang Observatory (BBO) which 
is so far motivated to detect the GW background produced during infiation. LISA, LIGO 
and BBO will be able to probe part of the history of the universe in the temperature range 
100 GeV-10'' GeV. The GW signal coming from particle physics phase transitions is directly 
related to the scalar potential describing the evolution of the order parameter. Observation or 
non-observation of GW will allow to put constraints on the parameters of these potentials. 
The measurement of the GW spectrum (peak frequency and intensity) can discriminate 
among different models (once combined with experimental measurement at colliders, for 
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instance, once knowing the higgs mass) and put constraints on the model parameters. For 
example, at LHC, we will be able to measure the Higgs mass but not the quartic or cubic self 
coupling of the Higgs. Only a linear collider can provide this information, which timescale 
could be beyond LISA. LISA could start constraining model parameters before a linear 
collider. In addition, LISA is sensitive to the 10 TeV scale which is beyond the reach of the 
near future collider experiments. 

The gravitational wave signal from phase transitions at around 10-100 TeV temperatures 
could entirely screen the signal from inflation, which detection is one of the main motivations 
for building BBO. 

We emphasize that our quantitative analysis can only be indicative given the uncertainties 
both at the experimental and theoretical level. The sensitivities of LISA, LIGO-III and BBO 
will certainly change during the next years. On the theoretical side, we use the estimate of the 
GW power spectrum of [6-9, 11, 12, 14] which is enough for the point of this paper. However, 
it certainly deserves improvement. We do not venture into this aspect in this work and just 
encourage that this question be re-examinated given the exciting experimental prospects for 
GW detection from phase transitions we demonstrated in this analysis. 

Note added 

As this work was being completed, Ref. [19] appeared where they re-examined the calculation 
of the gravitational wave background from turbulence. They disagree with and correct the 
dispersion relation used for gravitational waves in [11, 12, 14]. This leads, in particular, 
to a different prediction for the peak frequency as well as a different spectral dependence. 
Re-examination of the bubble collision spectrum is underway [20]. These new results can 
marginally affect the detectability regions of Figs. El ID and |S] but the overall conclusions will 
remain the same. 

Acknowledgments 

We are indebted to Alessandra Buonanno for very useful discussions. We also thank A. Nico- 
lis for a clarification on his work, M. Maggiore for pointing out some reference as well as 
C. Caprini and R. Durrer for discussions and finally James Wells for stimulating conversa- 
tions at the early stage of this project. 

References 



[1] URL http : //www . ligo . caltech . edu/ 



[2] URL http: //www. virgo . infn. it/ 



[3] URL http : //lisa ■ jpl . n asa . gov/| 



[4] T. L. Smith, M. Kamionkowski and A. Cooray, arXiv:astro-ph/0506422[ 
[5] URL universe .nasa.gov/program/bbo .html . 



16 



[6] A. Kosowsky, M. S. Turner and R. Watkins, Phys. Rev. Lett. 69, 2026 (1992). 

[7] A. Kosowsky, M. S. Turner and R. Watkins, Phys. Rev. D 45, 4514 (1992). 

A. Kosowsky and M. S. Turner, Phys. Rev. D 47, 4372 (1993) |arXiv: astro-ph / 92 1 1 004| . 

[9] M. Kamionkowski, A. Kosowsky and M. S. Turner, Phys. Rev. D 49, 2837 (1994) 
[arXiv:astro-ph/9310044| . 

[10] K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, Phys. Rev. Lett. 
77, 2887 (1996); K. Rummukainen, M. Tsypin, K. Kajantie, M. Laine and M. E. Sha- 
poshnikov, NucL Phys. B 532, 283 (1998); F. Csikor, Z. Fodor and J. Heitger, Phys. 
Rev. Lett. 82, 21 (1999). 

[11] A. Kosowsky, A. Mack and T. Kahniashvih, Phys. Rev. D 66, 024030 (2002) 
i arXiv: astro-ph /0 1 1 1483| . 

[12] A. D. Dolgov, D. Gr asso and A. Nicohs, Phys. Rev. D 66, 103505 (2002) 
|arXiv:a stro-ph /020646JJ . 

[13] R. Apreda, M. Magg iore, A. Nicohs and A. Riotto, NucL Phys. B 631, 342 (2002) 
|arXiv:gr-qc/0107033i . 

[14] A. Nicohs, Class. Quant. Grav. 21, L27 (2004) [arXiv:gr-qc/0303084l . 

[15] R. Bahantini et al, |arXiv:gr-qc/Q502054l 



[16] M. Maggiore, Phys. Rept. 331, 283 (2000) |arXiv:gr-qc/990900ll. 

[17] A. J. Farmer and E. S. Phinney, Mon. Not. Roy. Astron. Soc. 346, 1197 (2003) 
I arXiv:astro-ph/0304393 1 . 



[18] C. Cutler and J. Harms, Phys. Rev. D 73, 042001 (2006) |arXiv:gr-qc/0511092 



[19] C. Caprini and R. Durrer, |arXiv:astro-ph/0603476[ 
[20] C. Caprini, R. Durrer and G. Servant. In preparation. 
[21] L. Randall and G. Servant. To appear. 
[22] P. J. Steinhardt, Phys. Rev. D 25, 2074 (1982). 

[23] C. Grojean , G. Servant and J. D. Wells, Phys. Rev. D 71, 036001 (2005) 
jarXiv:hep-p h/0407019| . 

[24] A. Buonanno, G. Sigl, G. G. Raffelt, H. T. Janka and E. Muller, |arXiv:astro-ph/0412277| 



[25] C. Delaunay, C. Grojean, G. Servant and J. D. Wells, in progress. 
[26] J. R. Espinosa, M. Losada and A. Riotto, arXiv: hep-ph /04090 70l 



17 



